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Abstract 

\Q , 

^S) , We quantize subcritical bubbles which are formed in the weakly first order phase 

transition. We find that the typical size of the thermal fluctuation reduces in the 
quantum-statistical physics. We estimate the typical size and the amplitude of 
Q\ | thermal fluctuations near the critical temperature in the electroweak phase transi- 

tion using quantum statistical average. Furthermore based on our study, we give 
^2 ' implication on the dynamics of phase transition. 
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Ph! 1 Introduction 

a: 

The scenario of baryogenesis in the electroweak scale |J is attractive because we can get 
the data of this energy scale from recent experiments of elementary particle physics. How- 
ever, when one trys to explain the present baryon to entropy ratio in the minimal standard 
model, some serious problems exist. l)In order to generate sufficient baryon number, we 
need a deviation from thermal equilibrium. This can be solved if the phase transition 
is strongly first order with supercooling. However, in the minimal standard model, elec- 
troweak phase transition is estimated to be weakly first order in the 1-loop approximation 
@. 2) CP- violation source from CKM phase is too small to explain present baryon to 
entropy ratio. 3) We could not avoid washing out of baryon number by the sphaleron 
transition after the phase transition in the minimal standard model. Improvements and 
more precise analysis are now in progress to solve the above problems. In this paper, we 
devote to the analysis of the phase transition. To understand the dynamics of the elec- 
troweak phase transition, one may need to consider a behavior of the thermal fluctuation 
of the scalar field. We study on a thermal fluctuation using subcritical bubbles which was 
first suggested by Gleiser et al[|J. In these analysis bubbles are treated as the classical 
object and assumed that the shape is spherical. 

Recently one of the present authors insisted that the fluctuation of the subcritical 
bubble should be treated quantum mechanically. So one must quantize the bubble and 
consider up to higher energy levels. In this situation, one should take into account the 
modification from the 0(3)-symmetry of the bubble because the energy difference between 



each level is almost the same as differential energy. Another reason that we need to take 
into account the above modification exists. In our previous studies we have convinced 
that the phase mixing is sufficiently attained in the experimentally allowed region on the 
minimal standard model. From the above reasons, the modification cannot be neglected. 

This paper is organized as follows. In Sec. II we summarize our previous studies. 
In Sec. Ill we quantize the collective coordinate of the subcritical bubble including the 
deviation form 0(3)-symmetry and estimate the quantum-statistical averaged radius of 
the bubble. In Sec. IV we calculate the field amplitude inside the bubble. In the final 
section, we summarize our study. 



2 The typical size of classical subcritical bubbles 

We summarize in this section our previous results on the classical bubbles. In the weakly 
first order phase transition, the time development of the phase transition strongly depends 
on the thermal fluctuation. In the previous study M, we have assumed 



(/> + exp 



R?{t) 



for the spatial profile of the thermal fluctuation. + is the field value at the asymmetric 
vacuum phase. The Hamiltonian for the radius R(t) becomes 

H(p ' R) = mR) + l M(H) ' (2) 

where 

This profile will be justified near the critical temperature at which two vacua degenerate 
because almost all bubbles have thick walls. Here we neglected the term from the potential 
of the scalar field because that term vanishes near the critical temperature. Thus we 
obtained the typical size of thermal fluctuation by thermally averaging the radius; 

JdRdPRe W [-PH(P,R)} _ 2y^T 
{ )T fdRdPexp[-f3H(P,R)} ~ U 

Note that this is smaller than the correlation length. On the other hand, Gleiser et al. 
have assumed that the size of the bubble is equal to the correlation length[|J. We think 
that the typical size should be determined from the dynamics of the phase transition. 
However, we encounter with a serious problem: Can one treat the thermal fluctuation as 
the purely classical object? Actually, the number of state inside the bubble is too small; 

n(r < (R) T ) ~ 0(1) (5) 

Furthermore the condition for the decoherence, another measure for the system to be 
classical, is given by 



277T 3 / 2 \ 1/3 

(R)t <Rc -=[ ~r — n ~ 0.084GeV _1 . (6) 
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On the other hand, (R)t ~ 0.012GeV~ 1 , clearly, (R)t < R c - Thus, one cannot expect 
that the subcritical bubble is a classical object. 

From the above consideration one of the authors proposed that the fluctuation should 
be treated as quantum one with statistical fluctuation- dissipation in weakly first order 
phase transitions 0. In this paper, we first treat the subcritical bubble as purely quantum 
and then we statistically average the quantum bubbles. If we obtained, after our present 
analysis, no drastic deviation from the previous estimate, we would conclude again, with- 
out solving the Master equation for the density matrix, that the phase transition cannot 
accompany the supercooling. 



3 The typical size of quantum subcritical bubbles 
3.1 Subcritical bubble with modification 

We notice that the energy difference is the same order of magnitude as the modification 
energy from the spherical symmetry in the statistical averaging. Hence, in advance, we 
consider the contribution of the modification from the spherical configuration. We express 
the spatial profile of the thermal fluctuation as 



0(x) = 5^0/, m (r)¥i iTO (fi), 
where Y^ m (fi) is the spherical harmonics. Further, we assume the thick wall 

4>Lm(?) = 0+exp 



(7) 



R 2 (t) 



for each components. Here <fr + denotes the vacuum expectation value of the scalar field 
in the broken phase. The Hamiltonian becomes 



H(P,R) = H + J2(H e -H { 



(9) 



where 



W*) = ^ + ?m(*) + ^ + i)m(*). 



(10) 



3.2 The quantization and the measure problem 

Let us quantize the above system. Expressing the arbitrariness in the operator ordering 
by the parameter b, we write the quantum Hamiltonian as follows; 

where a := 2M{R)/R and c t := [(2/5) + (16/15)£{£ + l)](M/R). 

In the following discussion, we set b/a = N and we concentrate on only two cases 
N = 1 and N = 2. Here one should be careful on the measure in the latter case. The 
former case N = 1 corresponds to the operator ordering; 

Hi(P,R) = P^^P + CiR, (12) 
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where P := — id R . Thus, the integral measure of this case is constant. The latter case 
N = 2 corresponds to 

Ht(2p,z) = -\% + --; (13) 
where z := (R/R ) 3 ^ 2 , R := (R/M) 1 / 2 and z p := — id z . Hence the integral measure is 

poo poo 

/ dzoc dRR l/2 . (14) 

JO JO 

3.3 The wave function 

Let us construct the wave function of the subcritical bubble. The Schrodinger equation 
is given by 

i^e(R,t) = H e V e (R,t). (15) 
at 

For the stationary state = ipee~ lEet , Heipe = Eftyt ), the Schrodinger equation becomes 



-x>; + 

where we have defined 



+ l)x~ +d e [x- 



2N y 2N ' \ 2d 



Xi = fiX* (16) 



X t :=R-™ij) t {R) and x := a 1/2 R (17) 

and 

In order to solve this equation, we take the approximation, 

(19) 

and we neglect the first term of the potential. Then the Schrodinger equation reduces to 

- x '' +d >(*-0 x '=£, x >- < 2o > 

This approximation will be examined in the later discussion. 
Now, we define 

f- Pi) 



and then eq. (20) becomes 



w +ex '=w Xt - (22) 



As this equation is the same as that for the harmonic oscillator, the eigenfunctions and 
eigenvalues of this equation are given by 

AV(0°c(-l) n ^^, (23) 



and 



t Ln 



Ad 



3/2 



2n + l. 



(24) 



respectively. 

Using the above eigenvalues, we examine the approximation used on the potential 
under the assumption of fllTf). The original potential before the approximation is 



Ve(x) 



^(^ + l)^ 2 + ^-x ) 2 



(25) 



where Xq := et >n /2dg. Let us call the minimum value at which the first term coincides 
with the second term to be x*. 
In the case of N = 1 



1 + \ 1 + 



2 • 3V2 



for n = 0, 1. 



(26) 



2n + 1 

This is larger than xq and the approximation brakes down. On the other hand, for n > 2, 



i - wi 



2 • 3V2 
2n+ 1 



(27) 



is smaller than xq and it means that the approximation is reasonable. In the case of 
N = 2, in a similar way, we can show that the approximation is reasonable only for 
n > 1. 

Thus, for low energy states, we have to use different approximation. The simplest 
remedy for that is as follows. We first expand the potential 



V(x) ~ V{x) + -V"{x){x - x) 2 + 



where x is determined by V'(x) = 0. Defining 



1 



V"(x) 



1/4 



the Schrodinger equation for the stationary state becomes 



d 2 Xo 



di] 2 



+ V 2 X e 



1 



\V"{x) 



1/2 



it, ~ V ^ 



The solution is given by 



and 



n 2 d n e~ v 
X e , n ( V ) oc (-l) n e v - 



dr] n 



\V"{x) 



1/2 



4d f 



V(x) 



2n + 1. 



(2f 



(29) 



(30) 



(31) 



(32) 



Under the above two approximation, the wave function badly behaves at R ~ 0. This 
violates the hermiticity of the Hamiltonian. However, this bad behavior simply comes 
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from our approximation and should not be taken seriously. Actually we can easily see the 
true behavior near R ~ 0. At x <C 1, the eq. ( |T^ ) becomes 

d 2 X t 1/1 \ o , s 



cfe 2 2iV V2iV 
Assuming Xe oc x n , this becomes 



"- 1 "2ivJl" + 2ivJ= - < 34 > 

Thus, the regular solution is oc x 1+1 / 2N oc i? 1+1 / 2Jv . The original wave function is 
given by tpi(R) oc R 1+1 ^ N . The first derivative is ip'(R) oc R}/ N . Therefore actually the 
Hamiltonian is a hermitian operator. 

3.4 The typical size 

Now, we can calculate the typical size by the quantum-statistical averaging; 



where 



Q J °°dRR a R\^ n (R)\ 2 
{H)t >»- JS dRR"\^ n {R)\ 2 ' { ] 

R a being the integral measure; for N = 1, a = 0, and for N = 2, a = 1/2. The 
typical size of subcritical bubbles does depend on the integral measure, however after our 
approximation of the potential (0), the dependence on N disappears. Then, the effect 
of the integral measure only appears in the wave function with small n. Fortunately 
even for these wave functions, we can easily show that this effect to the typical size is 
small enough. We now estimate the typical size of a bubble using the effective potential 
in a finite temperature with the 1-loop approximation [0. We use mw = 80.6GeV, 
rriz = 91.2GeV and m t = 174GeV for the M^-boson, the Z-boson and the top-quark 
masses. And we set the mass of Higgs; m# = 65GeV that is the lower limit from the 
experiment 0. Accordingly, + = 49.4GeV and the critical temperature: T c = 98.3GeV. 
The result of our numerical calculation for the typical size is, 

(R) Q ~ O.tmGeV- 1 (37) 

near the critical temperature. This size is larger than that of thermal average: (R)t ~ 
0.012GeV _1 . 

To clarify the difference between classical and quantum regime, we calculate the dis- 
tribution function(W / (i? c )); 

/oo 
dPR«e- tRAP W(P, R c ) 
-oo 

/oo rco / | | \ 

dPR a c / dR A e- lR * p p R c + -R A , R c - -R A 
-oo J—oo V Z A J 

= R?p(R c ,R c ), (38) 
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where p is the density matrix given by 

P (Rc + \Ra, Rc - \Ra) = p(R, R') = Tr e , 

Under the present approximations, the final form of the distribution function does not 
depend on the measure; 

W(R C ) =R«p(R c ,R c ) = R c Tr e , 

We depicted profile of quantum distribution function W(R) (solid line) and the classical 
distribution function (broken line) in figure 1. The location of the peak of W(R) is larger 
than that of classical distribution function. This is because the fluctuations in quantum 
distribution function comes from quantum fluctuation as well as the statistical fluctuations 
which is expressed in the classical distribution function. 



reJR')^,n(R) 



(39) 



-^XfJR) 



(40) 



4 Fluctuation strength in fields 

While, fixing R to the above averaged value, we can estimate the fluctuation in the field 
value 0. Following the manner of previous studies [|3| @, we set 



A(t) exp 

In this case the Schrodinger equation becomes 

/ 1 d 2 



(R) 2 



QJ 



where 



m 



2m dA 2 
3 



+ f t A z 1> t = E t iPe, 



(R)l 

and the wave function becomes 
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+ 2£(£+l) 



ip e Ja) oc (-l) n e 



m 



,2 d n e 



V2 



da r > 



where we define 

a=(2mf e )*A. 
The fluctuation of can be calculated naively by 



where 



[<P 2 ) = 



(.a 



J dAA 2 ipe^ n (A) 



JdAij e , n (A) 

The result of our numerical calculation becomes 



30.4GeV. 



(41) 

(42) 
(43) 

(44) 
(45) 

(46) 

(47) 
(48) 



At the reflection point, = 0* = 10.5GeV. The fluctuation (|48| ) well exceeds 0* near 
the symmetric vacuum. This fact means that the true and false phases are well mixed at 
the scale typical in the fluctuations above the critical temperature. Therefore the use of 
1-loop approximation of the effective potential is no longer valid. 
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5 Summary 



In this paper, we quantized the subcritical bubbles and using quantum statistics, we 
calculated the typical size of subcritical bubbles as well as the typical fluctuation of fields. 

As a result of our calculations, we found that the typical size of subcritical bubbles 
is larger than that calculated by the classical thermal average, and that the size of the 
typical fluctuation of the field similarly well exceeds the first reflection point in our case. 
Therefore we may conclude that the phase transition cannot accompany supercooling. On 
the other hand, we must confess that the above study is quite qualitative. Furthermore, 
we could not justify whether the electroweak phase transition is first order or not as long 
as we use the 1-loop approximation for the effective potential, not effective action. 

Finally, we must mention the assumptions used here. First, we have assumed the 
Gaussian distribution for the spatial profile of the thermal fluctuation as in the previous 
studies. However there is no dynamical justification on this Gaussian distribution. Second, 
in order to consider the deviation from the spherical symmetric configuration, we assumed 
that each component does not depend on the parameter (£, n). We hope that we can justify 
the above assumptions in our quantitative study in the near future. 
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1. Fig.l: Quantum statistical distribution function as a function of the radius R of the 
subcritical bubble (solid line). Classical distribution function is also shown (broken 
line). The unit of the horizontal axis is GeV" 1 . 
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